Exploring the thermodynamics of a two-dimensional Bose gas 
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Using in situ measurements on a quasi two-dimensional, harmonically trapped 87 Rb gas, we infer various 
equations of state for the equivalent homogeneous fluid. From the dependence of the total atom number and the 
central density of our clouds with the chemical potential and temperature, we obtain the equations of state for 
the pressure and the phase-space density. Then using the approximate scale invariance of this two-dimensional 
system, we determine the entropy per particle. We measure values as low as 0.06 &b in the strongly degenerate 
regime, which shows that a 2D Bose gas can constitute an efficient coolant for other quantum fluids. We also 
explain how to disentangle the various contributions (kinetic, potential, interaction) to the energy of the trapped 
gas using a time-of-flight method, from which we infer the reduction of density fluctuations in a non fully 
coherent cloud. 
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Physical properties of homogeneous matter at thermal equi- 
librium are characterized by an equation of state (EoS), i.e. a 
relationship between some relevant state variables. For a fluid 
of particles, possible EoS's consist in expressions of pressure, 
density or entropy as functions of temperature T and chemical 
potential /i. For an ideal gas the EoS can be calculated exactly 
in any dimension for a Bose or Fermi gas. In the presence of 
interactions, one has to resort to approximations or numerical 
calculations, and a comparison with experiments is crucial to 
test their validity. Trapped atomic gases at thermal equilib- 
rium provide a powerful tool for this purpose [ 1 ] . Within the 
local density approximation, any intensive state variable takes 
at a point r in the trap the same value as in a homogenous 
system with the same temperature and the shifted chemical 
potential /i — V(r), where V(r) is the confining potential. 

The case of an interacting two-dimensional (2D) Bose fluid 
is particularly interesting in this context. Firstly, at non-zero 
temperature the Mermin-Wagner theorem precludes Bose- 
Einstein condensation I210. Therefore the EoS is expected 
to be continuous at any point, in spite of the existence of 
a superfluid, infinite-order phase transition, which is of the 
Berezinskii-Kosterlitz-Thouless (BKT) type [4 5|. Secondly 
the EoS of a 2D Bose fluid is scale invariant (6) in the so- 
called quasi-2D regime. The later refers to the experimen- 
tally relevant situation where single-particle motion is frozen 
along one axis, making it thermodynamically 2D, but where 
collisions still keep their 3D character Q. The (approximate) 
scale invariance originates from the fact that in this regime, 
the coupling strength g is an energy-independent dimension- 
less coefficient, and thus provides no energy, nor length scale, 
in contrast with the ID or 3D cases. This implies in particular 
that dimensionless thermodynamic variables such as the phase 
space density T> or the entropy per particle S are functions of 
the ratio n/k^T only, with g as a parameter. 

Recent experiments with trapped 2D Bose atomic gases 
have demonstrated the existence of a BKT-type transition at 
a threshold phase space density. One line of investigation ex- 
ploited matter-wave interference to monitor the appearance of 
an extended coherence in the sample |8] [5), and another ap- 



proach used a time-of-flight (ToF) technique to measure the 
momentum distribution of the gas ifTOll . The steady-state scale 
invariance was verified in ifTTI . In this paper we present a 
detailed experimental investigation of several thermodynamic 
properties of a 2D Bose gas. We describe measurements of 
the EoS for the pressure from a count of the total atom num- 
ber of a trapped 2D Bose gas, with a wide range of thermo- 
dynamic parameters. From the same set of data we use the 
central spatial density to access the EoS for phase space den- 
sity. Combining these two measured EoS with the scale in- 
variance we obtain the EoS for the entropy per particle. We 
show that this quantity rapidly decreases around the super- 
fluid transition point and then approaches zero in the highly 
degenerate regime. We also present an original method to ex- 
tract from a ToF in only one direction of space, the various 
contributions (kinetic, potential, interaction) to the total en- 
ergy of the trapped gas. This method is applicable to any low- 
dimensional fluid. Here it shows that density fluctuations of 
our 2D Bose gas are essentially frozen even when its thermal, 
non coherent fraction is significant. 

Our 2D Bose gases are prepared along the lines detailed in 
lfT2ll . We start with a 3D Bose-Einstein condensate of 87 Rb 
atoms confined in a magnetic trap in their F = mp = 2 in- 
ternal ground state with an adjustable temperature. We slice a 
horizontal sheet of atoms with an off-resonant, blue-detuned 
laser beam with an intensity node in the plane z = 0. It pro- 
vides a strong confinement along the direction perpendicular 
to this plane, with oscillation frequency u> z /2tt = 1.9 (2) kHz, 
which correspond to the interaction strength g = \^8ira/£ z f=s 
0.1, where a is the 3D scattering length and £ z = y/h/mu) z . 
The energy fiuj z is similar to or larger than the thermal energy 
k^T and the interaction energy per particle, so that most of the 
atoms occupy the ground state of the vibrational motion along 
z. The magnetic trap provides a quasi-isotropic confinement 
in the xy plane with frequency lu j 2tt — 20.6 (1) Hz IT31 . 

After an equilibration time of 3 seconds in the combined 
magnetic+laser trap, we measure the in situ density distri- 
bution of the gas by performing absorption imaging with a 
probe beam propagating along the vertical axis. The conven- 
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FIG. 1: (Color on line) Absorption imaging of quasi-2D clouds of 
87 Rb atoms, (a) Image obtained with a short pulse (~ 2/xs) of an 
intense probe beam (I/I sat = 40). (b) Image obtained with a longer 
pulse (50 fis) of a weak probe beam (J// sa t = 0.5). (c) and (d) 
Radial density profiles for image (a) (hollow circles o) and image (b) 
(filled circles •) in linear (c) and logarithmic (d) scales. 



tional procedure where one uses a weak probe beam with an 
intensity I well below the saturation intensity 7 sat is prob- 
lematic in this context lfT2l . Indeed for the relevant range 
of temperatures (40-150 nK), the atomic thermal wavelength 
At is comparable to the optical wavelength used for probing, 
A op t = 780 nm. Consequently in the highly degenerate re- 
gion of the gas (V = n' 2D ^A T ^> 1), the average distance 
between neighboring atoms is much smaller than A op t and the 
absorption of a weak probe is strongly perturbed by collective 
effects. To circumvent this problem we probe the gas with a 
short pulse (duration ~ 2 /is) of an intense probe beam (typi- 
cally ///sat = 40 to 100) 03). The interaction of any given 
atom with light is then nearly independent of its neighbors. 

High-intensity imaging, which was also used in [11 1, pro- 
vides a faithful measurement of the atomic distribution in the 
central region of the trap, where the density is large. How- 
ever the quality of the images suffers from a large photon shot 
noise, which spoils the detection of the low-density regions 
of the cloud (Fig. [T^). In order to probe reliably these re- 
gions on which we base our determination of temperature and 
chemical potential, we complement the high-intensity imag- 
ing procedure by the conventional low -intensity one (Fig.[lJ)). 
In practice for any set of parameters to be studied, we perform 
one run of the experiment with high-intensity imaging, and 
one with low-intensity imaging immediately after. The repro- 
ductibility of the experiment is checked by acquiring several 
pairs of images for a given set of experimental parameters. 

The procedure for image processing is detailed in the Aux- 
iliary Material. In short, for each pair of images it provides 



the temperature T, the chemical potential at center /i and the 
density n(r) at any pixel of the image. We assume the atoms 
in the excited states of the z motion to be described by the 
Hartree-Fock mean-field (HFMF) theory ifTOl [B1JT7) : there- 
fore, we can self-consistently calculate the populations of the 
excited states, and subtract it from n(r) in order to obtain the 
density distribution i%o(r) in the ground state. The validity of 
this procedure was checked by analyzing the results of a quan- 
tum Monte Carlo calculation for a range of parameters similar 
to ours [ 18 1. 

We start our thermodynamic analysis by inferring the pres- 
sure T) of the homogeneous gas from our measure- 
ments. Here we adapt to the two-dimensional case the tech- 
nique presented in JT] and that has successfully been used 
in 3D for Fermi gases |fl9) . We show that P(fi,T) is di- 
rectly related to the atom number N$ = J no(r) cPr in our 
harmonic trap. Indeed, the local density approximation re- 
lates the density n$(r) to the density of the homogenous 
gas [m — ^( r )j T]. For an isotropic harmonic potential 
V(r) = muj 2 r 2 /2 the total atom number is 

IULU J_ oc 

and using the thermodynamic relation n[^l — (dP/dfi) T , 
we find Ao = (27r/mcj 2 )P(/i, T). Introducing the dimen- 
sionless quantity V = P\\/k-QT, which we refer to as the 
reduced pressure, we then obtain 

where u is to be replaced by the geometrical mean of u x and 
uj y for an non-isotropic potential. Our results for the pressure 
are summarized in Fig. [2^, where we plot V deduced from 
Eq. (j2j) as function of ^jb/k^T. The temperatures of the data 
entering in this plot range from 40 nK to 150 nK. The fact 
that all data points collapse on the same line show that V is a 
function of the ratio /i/k^T only, as expected from the scale 
invariance of the system. The HFMF theory is represented by 
a continuous line in the normal region and by a dotted line in 
the superfluid region. The dashed line is the Thomas-Fermi 
prediction at zero temperature V = TT(fi/k^T) 2 /g. The grey 
area is the parameter subspace accessible for an ideal Bose 
gas. Interestingly, although the phase space density V can 
take arbitrarily large values in an ideal 2D Bose gas, one can 
show that the reduced pressure is < 7r 2 /6, where the equality 
provides a local criterion for Bose-Einstein condensation in a 
trapped ideal gas. 

We show in Fig. our measurements for the phase space 
density T>, obtained from the central density of each cloud. In 
wide gray line we plot the prediction of |6|, which is in good 
agreement with our results. A measurement of T> was also 
reported in [11] for a quasi-2D Cesium gas, showing a simi- 
lar agreement with J6). To be quantitative we have fitted the 
prediction of [6 1 to our data by multiplying it by a global fac- 
tor £, and obtained £ = 0.93 as optimal parameter. This 7% 
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FIG. 2: (Color on line) Equations of state for (a) the reduced pressure V, (b) the phase space density T> and (c) the entropy per particle S. The 
Hartree-Fock mean field prediction is plotted in full line and extended in dotted line beyond the expected superfluid transition. The dashed line 
indicates the Thomas-Fermi prediction. In (a) the grey area indicates the region of parameter space accessible to an ideal gas. In (b) the thick 
grey line indicates the prediction from [6|. 



discrepancy may be due to residual loss of detectivity in high 
density regions. Note that the measurement of the pressure 
EoS is much less sensitive to this possible bias since it relies 
on the count of the atom number iVo over the whole cloud 
rather than on the highest value of the spatial density. 

From our measurements of V and T> we also obtain the 
equation of state for the entropy per particle <S(/i, T): 

- = 2^-^, (3) 
k B T) k B T ' 

which can be derived starting from the entropy per unit area 
s = (dP/dT) , assuming the EoS for V to be scale invari- 
ant l20l . The corresponding result is shown in Fig. [2};. As 
expected, S is large in the non-degenerate regime and rapidly 
decreases around [ijk B T w 0.15, where the superfluid tran- 
sition is expected for our value of g ETj . Finally S tends 
to zero in the Thomas-Fermi regime. Our data points with 
the largest phase-space density (fj,/k B T > 0.5) correspond to 
S = 0.06 (1) fee only. Note that since the BKT transition is of 
infinite order, one does not expect any discontinuous change 
for V, T> or S at the superfluid transition for an infinite homo- 
geneous fluid, although the superfluid density jumps suddenly 
from0to4/A| E2. 

We now turn to the last part of our study, where we il- 
lustrate how to measure the various contributions to the en- 
ergy of our trapped 2D gases: potential energy E p in the ex- 
ternal trapping potential, kinetic energy of the particles E^, 
and interaction energy between atoms E\. We first point out 
the simple relation E p = E^ + E- l7 obtained from virial 
theorem assuming 2D contact interaction. We can measure 
E p = J no(r) V(r) d 2 r, from an in situ image, but we still 
need to disentangle the contributions of E^ and E\ to the total 
energy. This can be done by abruptly switching off interac- 
tions at time t = 0. Each particle then undergoes a free har- 
monic motion r(t) — cos(ut) r(0) + sin(u;t) v(0)/ui. The 
potential energy after a time t following the switching off of 



the trapping laser is given by 

E p (t) = E p (0) cos 2 (wt) + E k {0) sin 2 (ujt), (4) 

where we used the fact that the correlation (r(0) • v(0)) is 
zero at thermal equilibrium. Thus we can extract £"k(0) from 
the time evolution of E p , which we obtain from the density 
profiles at different times t. 

In order to implement this procedure, we perform a "one- 
dimensional" ToF by switching off abruptly the laser provid- 
ing the confinement along z while keeping the magnetic con- 
finement in the xy plane. The gas then expands very fast along 
the initially strongly confined direction z, as shown in figures 
[3^ to |3p, and interactions between particles drop to a negli- 
gible value after a time of a few lu^ 1 , where wj 1 ^100 fxs. 
The subsequent evolution in the xy plane occurs on a longer 
time scale given by a; -1 ~ 8 ms. From Eq. Q and £k(0) < 
E p (0), we expect the size of the gas to decrease for t < a; -1 , 
which can be understood in simple physical terms. The equi- 
librium state of the 2D gas results from a balance between 
the trapping potential, which tends to compress the gas, and 
the kinetic and interaction energies, which tend to increase its 
area. When interaction energy drops to zero the equilibrium 
is broken and the gas implodes in the xy plane. A similar ID 
ToF technique was used recently in Boulder with the value 
of t fixed at 7r/2w iflOl . For this particular choice the initial 
momentum distribution is converted into position distribution 
and can thus be measured accurately ll23l . 

We show in figure^ an example of measurement of E p (t) 
for a gas with N = 6.1 10 4 , T = 72 nK, and ^/k B T = 0.59. 
From the contraction of the gas, we infer E^/E p — 0.56 (3), 
from which we deduce E-JE P — 0.44 (3) using virial theo- 
rem. This configuration is thus neither completely in the very 
dilute regime (E\ <C E\<_ ~ E p ) nor in the Thomas-Fermi 
regime (E k <C E\ ~ E p ) and contains comparable thermal 
and quasi-coherent fractions. 

The measurement of E\ is of particular interest in this 
case since it gives access to the density fluctuations in the 
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FIG. 3: (Color on line) (a) to (d) Side view of a cloud initially in 
the 2D regime and expanding along z once the laser providing the 
confinement in this direction has been switched off. (a) t — 1 ms; 
(b) t — 2 ms; (c) t = 3 ms; (d) t = 4 ms. (e) Time evolution of the 
potential energy E p . The different lines represent a fit to the data of 
a parabola (solid black line), the time evolution assuming completely 
frozen fluctuations (dashed red line) and the one expected for a dilute 
non-condensed gas (dash-dotted green line). 



gas. Indeed, by definition E- t = (h 2 g/2m) §{n^{r)) d 2 r = 
{fi 2 g/2m)J r J (no(r)) 2 d 2 r, where we have introduced the 
parameter T which characterizes the degree to which density 
fluctuation are reduced. In the limiting case of a very dilute, 
non-condensed gas, one expects T = 2, since (tiq) = 2(no) 2 , 
while in the opposite limit of suppressed density fluctuations 
J- = l. Since our measurement provides us with E- u we can 
infer the value of T, from the comparison with the quantity 
(h 2 g/2m) J (no(r)) 2 d 2 r, calculated using the in situ density 
profile no. For the experimental conditions of figure [3£, we 
find T = 1.1 (1), very close to 1 that would correspond to 
completely frozen density fluctuations. Note that this is ob- 
tained for a gas still far from the Thomas-Fermi limit since 
i?k ~ El This "early" freezing of density fluctuations is 
an important ingredient for the proper operation of the BKT 
mechanism. This presuperfluid phase, whose existence was 
also inferred by different methods in ifTUl and ifTTl . consti- 
tutes a medium that can support vortices, which pair at the 
superfluid threshold. 

In conclusion we have presented in this Letter various as- 
pects of the thermodynamics of a 2D Bose gas, investigating 
first the EoS's for the pressure, the phase space density and 
the entropy. Our results confirm the scale invariance that was 
discussed theoretically in |6 1 and observed in [ 1 1 1 for T>. We 
point out that the entropy per particle drops notably below 
0.1 fee beyond the transition point. With such a low entropy 
these 2D Bose gases can constitue excellent coolants for other 
quantum fluids such as 2D Fermi gases [24 J. We have also 
presented a method that allows one to extract the various con- 
tributions to the total energy of the system. By applying it to a 
degenerate but not fully coherent 2D cloud, we find that den- 



sity fluctuation are nearly frozen, marking the presuperfluid 
phase. 
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Auxiliary material 

Imaging a dense 2D atomic cloud. The calibration of ab- 
sorption imaging consists in relating the number of missing 
photons on a pixel to the number of atoms on this pixel. The 
interaction between a probe beam and a single atom is char- 
acterized by the absorption cross section a defined by the re- 
lation 7 = crl /(Tuo-l), where 7 is the photon scattering rate, 
/ the intensity of the beam on the atoms and wl/27t its fre- 
quency. In the case of a monochromatic resonant beam prob- 
ing a two-level atom 



I 



7 



2 1 + L, 



(5) 



In the limit where / <C / S at the absorption cross section is 
cr = Thoj L /2I sat . 

In practice one must take into account stray magnetic fields, 
non-zero linewidth of the probe laser, optical pumping effects, 
etc. To model this complex situation, we heuristically replace 
I sat by an effective saturation intensity al sa t and T by an ef- 
fective linewidth (3T. We then write the number of photons 
iV p scattered during an imaging pulse of given duration r 



AT - ? V 1 

iV n = 7T — T, 

p ' 2 I + aI sat 



or equivalently 



a = cr - 



(6) 



(7) 



At low intensity N p is proportional to I as in the two-level 
case, but with a multiplicative coefficient /3/a due (for exam- 
ple) to the broadening of the resonance line. At large intensity 
the number of scattered photons saturates at (3Tt/2 instead of 
r-r/2, which models a reduction that can be caused by optical 
pumping effects, for instance. 

We now turn to the description of absorption imaging of a 
2D atomic cloud. The imaging process consists in shining a 
resonant laser beam on an atomic sample, and in imaging the 
transmission of the sample on a camera. In order to relate the 
missing photon number to the atomic density n(x, y), we cal- 
culate the probability for a photon of the probe beam to reach 
a pixel of the camera. We introduce the area A associated to 
this pixel in the atomic plane. In the limit where a <C A a 
photon has a probability a /A to be absorbed by a given atom, 
hence aprobability P t = (l — a/A) N to be transmitted, where 
N is the number of atoms in the area A. Thus we find: 



Pt 



(8) 



where we have used n — N/A, assuming that the atomic den- 
sity varies smoothly over the pixel size. The intensity of the 
beam at the output of the cloud is If = P t Ii and we obtain: 



In 



I({x,y) 
h 



(9) 



where a depends on the effective intensity / on the atoms [Eq. 
Q]. If the optical thickness of the cloud is large, i.e. if the 
intensity If just after the plane of atoms is significantly lower 
than the intensity just before this plane, the effective in- 
tensity I must be determined in a self-consistent manner by 
imposing: 



If = Ii-n a(I) I. 



(10) 



The elimination of the effective intensity / from Eqs. (f7|)-(|T0]> 
yields: 



no~Q j3 



h - If 



(11) 



It is interesting to note that even though the derivation in a 2D 
system differs from the 3D case, the result is similar to the one 
given in fl4l . The first member of the right-hand side of Eq. 



( 1 1 1 is dominant in the weak intensity limit, and corresponds 
to the 2D analog of the 3D Beer-Lambert law. In the high in- 
tensity limit, the second member of the right-hand side dom- 
inates. We calibrated a = 2.6 (3) using the same method as 
in 1 14 1 : we performed absorption imaging of clouds obtained 
in similar experimental conditions with various intensities I{ 
ranging from 0.1 / sa t to 6/ sa t, and imposed that these mea- 
surements provide the same result for the left-hand side of Eq. 
( [TT| ). Here we restricted ourselves to low atomic density re- 
gions, to ensure that collective effects in the optical response 
of the gas were negligible. The calibration of (3 = 0.40 (2) 
was performed as in [ 12], using the HFMF prediction as a fit 
to the low-density parts of our atomic distributions, and using 
/1, T and /3 as optimization parameters. In [12] where only 
low intensity imaging was used, this calibration provided the 
detectivity factor £, which is related to the present parameters 
a and /3 by £ = (15/7)/3/a. 

Obtaining a density profile from an absorption image. 
The confinement potential in the xy plane is essentially pro- 
vided by our magnetic trap, but it may also be affected by 
some imperfections in the intensity profile of the beam that 
freezes the z degree of freedom. These imperfections are re- 
vealed by looking at the center of mass oscillations x cm (t) 
and y C m(t), shown in Fig.|4^. Whereas the oscillation along 
the direction of propagation of the "freezing laser" (x) shows 
no deviation with respect to harmonic motion, the oscilla- 
tion along y is damped. This is likely caused by irregular- 
ities of the transverse intensity profile of the freezing laser. 
In order to cope with these defects we have abandoned the 
standard technique consisting in making angular average of 
the images to produce radial density profiles. Instead we 
take advantage of the separability of the potential in the xy 
plane: V(x,y) — mu!^.x 2 /2 + U(y), where U(y) accounts 
for the magnetic trapping potential and the irregularities of 
the freezing laser. We consider cuts of the measured density 
profile along the x direction, measured for various y^'s with 
i = 1, . . . ,q. In practice, we consider the q — 31 central 
lines of our images. We expect that two cuts corresponding 
to yi and j/2 coincide, provided we shift the second one by 
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FIG. 4: (Color on line) (a) and (b) Center of mass oscillations (hollow 
circles o) along x (a) and y (b). The red lines correspond to a fit with 
a sine (a) and a damped sine (b). (c) Reconstructed potential along 
the y axis (filled circles •) and a harmonic fit (red line). 

making the substitution moj 2 x 2 — > mtu 2 x 2 + U(y2) — U(yi). 
In practice we perform a least-square fit to optimize the su- 
perposition of the various cuts, taking the numbers U{yi) as 
parameters. We use a single set of U(yj) to fit a whole series 
of images taken at a given temperature. The robustness of the 
procedure is excellent, as shown in Fig. [4}}, where we give the 
reconstructed potential U(y), with bars corresponding to the 
statistical errors of the U(yjYs for various series of images 
acquired at different temperatures. 

Analyzing a density profile. For each configuration of a 
2D cloud, we first take a low-intensity image and then a high- 
intensity image in the following run. From the density profile 
of the low-intensity image, we determine the temperature T 
and the chemical potential /i by fitting the low density region 
with the HFMF prediction. Our fitting function takes into ac- 
count the residual excitation of the z degree of freedom. The 
high-intensity image provides the density profile n(r). 

Once T and /i are known, we self-consistently determine 
the population of the excited states using the method described 
in lfl0l[T31 . assuming the atoms in the excited states j > 1 of 
the z motion to be in the HFMF regime. In practice we restrict 
the analysis to the first ten levels. In order to give an estima- 
tion of the contribution of the various levels j > 1 to the 
total density, we show in Fig. [5^ numerical results obtained 
by applying this procedure to a numerically generated pro- 
file, produced using the prediction |6| with T = 100 nK and 
n/k-gT — 0.45. This temperature is on the high side of our 
experimental range, where the influence of the atoms in the 
excited states along z is expected to be the most important. 
We plot in Fig. [5^ the phase space density of the excited states 
2?( cxc ) 5 distinguishing the contribution of the state(s) j = 1, 
j = (1, 2), j = (1, 2, 3), etc. For comparison we also plot the 
profile 2?(°) obtained from (6), associated to the atoms in the 
ground state. Note that the contribution of the states j > 4 



is already negligible. The phase space density associated to 
each excited state is lower than 0.5, which justifies to treat the 
atoms in these states within the HFMF approximation. The 
flattened shape of the density distributions in the central re- 
gion is due to the repulsive interaction with the atoms in the 
ground state of the z motion. 

This procedure also allows us to calculate the effective po- 
tential felt by the atoms in j ~ 0, when the repulsive po- 
tential W(r) created by the atoms in j > 1 is taken into 
account. Plotting together W(r) and the trapping potential 
V(r) (Fig. |5|)) we see that W{r) is essentially negligible 
(< InK) and one can thus consider the density rto(r) to be 
insensitive to the presence of the atoms in j > 1. 

Effect of the finite interaction energy. The expression g = 
V&na/iz for the interaction strength assumes that the atoms 
are confined in the gaussian single-particle ground state of the 
harmonic motion along the z direction. However, because the 
interaction energy is not completely negligible compared to 
huj z , the state of the z motion at T = is modified by these 
interactions, which in turn modifies the coupling strength g. 
To estimate the corresponding effect we used first-order per- 
turbation theory to determine the modified ground state ifo(z) 
of the z motion. We then calculated the coupling strength g, 
which is proportional to J |yo| 4 dz. We find that g is reduced 
by the factor «1- 1.4n( 2D W 2 , which is rj 0.9 at the center 
of our densest clouds, comparable to the noise level on our 
data. 
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FIG. 5: (Color on line) (a) Phase space density of the ground (solid 
red line) and excited state(s) of the z-motion. The n-th line from 
the bottom corresponds to the contributions of excited levels 1 to 
n. (b) Comparison of the trapping potential (red solid line) and the 
repulsive potential created by the excited atoms on the population in 
the ground state (blue dotted line). 



